Fourier Analysis of Ghost Imaging 
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Fourier analysis of ghost imaging (FAGI) is proposed in this paper to analyze the properties of 
ghost imaging with thermal light sources. This new theory is compatible with the general correlation 
theory of intensity fluctuation and could explain some amazed phenomena. Furthermore we design 
a series of experiments to testify the new theory and investigate the inherent properties of ghost 
imaging. 
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Conventionally, coherent, at least partially coherent, il- 
lumination is required to image the phase detail of an ob- 
ject [1] , and the theoretical limit of the resolution is deter- 
mined by the wavelength A of the light beam. Since the 
coherent sources in the ultra-violation and even shorter 
frequency regions are still not existed, we could not ob- 
tain a higher resolution by using a coherent light with 
a shorter wavelength. In 1994, Belinsky and Klyshko 
[2] found that, by exploiting the spatial correlation be- 
tween two entangled photons generated by parametric 
down conversion (PDC), ghost imaging could be realized 
with entangled incoherent light. Due to the possibility to 
improve the resolution, ghost imaging is bestowed with 
great potential in quantum metrology, quantum hologra- 
phy and quantum lithography, etc. 

The nature of ghost imaging leads to many interesting 
debates about the necessity of the entanglement [3-14]. 
Presently, the idea that both classical thermal sources 
and quantum entangled beams can be used for ghost 
imaging and ghost diffraction is widely accepted. 

In order to realize practical use of the ghost imaging, 
further investigation about the properties of ghost imag- 
ing is stringent. Gatti et al. compared the visibility in 
the classical and quantum regimes, and pointed that the 
visibility of classical regime increases as the mean pho- 
ton number of per mode in thermal state increases [9] . 
D'Angelo et al. compared the resolution of quantum and 
classical ghost imaging [15]. Meanwhile, Ferri et al. real- 
ized the high-solution of ghost imaging and ghost diffrac- 
tion with thermal light [16], and they even tried ghost 
imaging with homodyne detection [17]. Bache et al. an- 
alyzed the spatial average technique which would make 
the imaging bandwidth of the reconstructed diffraction 
pattern virtually infinite and the correlation convergence 
rate faster [18]. Cheng et al. have presented theoreti- 
cal analysis on the noise of ghost imaging with entan- 
gled photons [19] and pseudo-thermal light fields [20]. 
Most recently our group demonstrated lensless Fourier- 
transform ghost imaging with both amplitude only and 



phase only object by using classical incoherent light [21]. 

There are still some questions unresolved, such as what 
limits the complexity of the object, why the visibility 
decreases as the transmission increases. In this paper we 
propose a new theory, Fourier analysis of ghost imaging 
(FAGI), to answer these questions. 

For an ideal infinite and uniform thermal source lo- 
cated at z = 0, the field at a point in the source is 
presented by U {x , y ,i) — J j ±oo dk x dk y Aexp[i(k x x + 
k y y — ujt)] . When k x + k y > (^) 2 ,fc z is imaginary, 
and the amplitude of this evanescent component decays 
exponentially as the distance z increases, so the inte- 
gral region can be confined in ki £ (—k,k),i e {x,y). 
Since the sources used in experiments are always fi- 
nite, which means k x and k y are not uniformly dis- 
tributed, the typical Gaussian distribution is adopted. 
For simplicity k x is assumed to be independent of k y , 
and the amplitude A is a time variant and random 
distributed, so the field is written as Uo(xo,yo,t) 

I f ±Qo dk x dk v exp(- fcx a t " 2 fca )Ag(t) exp^fc^o + k y y )}. 

Now we analyze the paradigmatic example of the imag- 
ing system given in Ref . [22] , except that the source is re- 
placed by a thermal source. Here the setup of the object 
arm is an f-f system with the object close to the beam 
splitter (BS). For simplicity the output-plane of the BS 
is selected as the source, the fields on these two planes 
are shown as 



U{x ,y Q ,t) 



dk x dk y exp( 



k'x + k y N 



±oo 



A^{t) exp[i(k x x + k y yo)} 



(1) 



V(x ,y 0) t) = - 



dk x dk y exp( o - / 



A%(t) exp[i(k x x + k y y )} 



(2) 
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and the object transmitting function is t(xo,yo). Ne- 
glecting the finite pupil of the lens, through the Fourier 
Transform of the lens, the spatial frequency distribution 
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at the focal plane of the object arm is 

U(x f ,y f ,i) = 



±00 



k 2 + k 2 
dk x dk y exp( — ^ — -) 



A s (t - -±)T[{—x f - k x ), (—y f - k y )} 

(3) 

Here, the constant phase, which has no sensible influ- 
ence on the result, is neglected. Because of the incoher- 
ence of the thermal source, a fixed pixel detector collects 
a certain range of spatial frequency components. Com- 
paring with a single frequency component in the coherent 
illumination case, it is easy to understand why the ghost 
imaging can be realized with only a pixel detector in the 
object arm. 

In this system, the reference arm is a 2f-2f setup, and 
the transmitting function h(x 2 f,y2f,x ,y ) = S(x 2 f — 



Xo,y2f — Vo), so the field distribution on the detector 
plane is 



k 2 x + ky 



V(x2f,V2f,t) = \ I I dk x dk y exp(- 

2 J J ±OQ a 



±00 

4/, 



A & - —) exv[i{k x X2f + k y y 2 f)] 



(4) 



The intensity fluctuation correlation function is 
G(xi,yi,X2,V2) = \T(xi, 2/1,^2,2/2) | 2 , and the sec- 
ond order correlation function T(xf,yf,X2f,y2f) = 
(U(xf,yf,t)V*(x2f,y2f,t)). On the assumption of the 
independence of the amplitude of different wave vectors, 
(A s (t - ^)A % , (t - ^)) = 0, when k ^ k' , the second 
order correlation function is simplified as 



J 



T(x f ,y f ,x 2 f,y2f) 



dk x dk y {A^(t - —)A % (t - -^-)) cxp(~2 kx+ 2 kv ) 

±00 C C (7 

2ir 2tt 



r 



(5) 



Within the coherent time, {AAt — —)Ar,(t — —)) is a const. When Xf = = ,eq.(5) is presented as: 



T(0,0 7 x 2f ,y 2f ) = -I 



k 2 + k 2 

dk x dk y exp(-2 x y )T(-k x , -k y ) exp[-i(k x X2f + k y y 2 f) 



(6) 



±00 



r 



this equation is a Fourier Transform, where the object properties of Fourier Transform, a convolution of the ob- 
spectrum multiplies a Gaussian function, which causes ject image and a Gaussian function is obtained: 
a increasing loss for high frequencies. According to the 



a I 

F(0, 0, x 2f ,y 2f ) = —t(-x 2 f, -V2f) ® exp(-cr 
lb 



2 x 2f + y 2f , 



(7) 



r 



The presentation of the Gaussian spreading function con- 
tributes a background to the object image, which ex- 
plains the occurrence of the inherent noise in the second 
order intensity fluctuation correlation. From this equa- 
tion it easy to deduce that when the transmission area 
increases, the background increases while the visibility 
decreases. 



When Xf^0,yj^0, the second order correlation 
function is 

2ir 

T(x f ,y f ,x 2 f,y2f) oc exp[-i — (x f x 2f +y f y 2f )}T(0,0,x 2f ,y 2f ) 

(8) 

If the pixel detector is replaced by a bucket detector 
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in the object arm, the corresponding result will be 



dxfdyfT(xf,Xf,x 2 f,y2f) oc J J dx f dy f exp[-i^(x f x 2 f + 2//2/2/)]r(0, 0, x 2 f, y 2 f) 



(9) 
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FIG. 2: The frequency response curve of Table. I. 



FIG. 1: (a)the image of different objects and (b)the corre- 
sponding spatial frequency spectrum, (i)two slits, the width 
is 300 ymi, and the distance is 900 yon. (ii)four slits, width 
150 fim, and distance 450 fim, (iii)six slits, width 100 fim, 
and distance 300 fxm. 



Similar results can be obtained with only coefficient 
difference. 

We have designed several experiments to demonstrate 
these deductions. Since the source we used is a pseudo 
thermal source, the best resolution is limited by the co- 
herence length l c , simply by replacing A with l c , all re- 
sults obtained in previous discussions are still correct. 
In the experiment, the coherence length of the source is 
75/ito . 

The first experiment is to investigate the frequency 
response of ghost imaging. The objects have the same 
duty cycle but different periods, and we get the image 
and the corresponding spatial frequency distribution of 
each object. 

In Fig.l, the black bold line in each (a) is the object, 
the disperse squares are the normalized intensity fluctua- 
tion correlation, and the fitting curve is the image. Each 
(b) is the corresponding spatial frequency distribution of 

(a) ' 

In this experiment we choose the ratio of the slit width 
and the distance a : d — 1 : 3 to make sure the image re- 



TABLE I: the relative frequency response rate in ghost imag- 
ing, (the coherence length is 75 fim) 



CSL a 



FFC 



RFR C 



infinity 
300 
150 
100 

75 





0.0625 
0.125 
0.1875 
0.25 



1 

0.64 
0.449 
0.388 




"Character Spatial Length. 
^Frequency of First-order Component. 
c Rclative Frequency Response. 



constructed by low frequency has no direct current back- 
ground. Moreover, comparing the first-order component 
of the image to the zero-order component of the object in 
Fig.l (i), the ratio is bigger than 1/2, which implies no di- 
rect current background of zero-order component, so the 
background noise is the contribution of the convolution 
of Gaussian spread function. As the frequency increases, 
the ratio of image first-order component to object first- 
order component decreases, so the effective information 
decreases. Results are shown in Tabled and Fig. 2. 

In Fig. 2 the frequency response curve has an approx- 
imate Gaussian form, as given in the theoretical result 
Eq.(6). 

In the following experiment we demonstrate the rela- 
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the transmission area is large enough, the image is sub- 
merged in the inherent noise. Here we design a two-slit 
and its reverse, the objects are shown in Fig. 4. 

in experiment, the object (a) is easy to be imaged, 
but under the same condition, its reverse (b) can not be 
imaged. 

From the above experiment results, we conclude that 
the Fourier analysis theory presented in this paper can 
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FIG. 3: The experiment result of objects with different trans- 
mission area, (a) the image of two-slit, (b) four-slit and (c) 
six- slit. 

tionship between the visibility and the transmission area. 
Here the objects have the same period parameters but 
different number of period, which suggests the spatial 
spectrums of each object have similar distribution and 
response. In eq.(7), the transmission area is composed 
by many points, this area can be expressed by a comb 
function, each point convoluted with a Gaussian spread 
function, and the image is the summation of all Gaussian 
spreading functions. As the transmission area increases, 
more points contribute to the background that directly 
makes the visibility decrease. The experiment result is 
shown in Fig. 3. 

The period d is 600/xto, and the width a of the slit is 
200/jto. From this experiment it seems reasonable to say 
the visibility decreases with the increase of the transmis- 
sion area. 

Based on the second experiment we suppose that when 



FIG. 4: (a)the object and (b)its reverse. 



interpret the principle of ghost imaging in different sys- 
tems, and resolve some difficulties of intensity correlation 
function. 

In summary, we propose a new theory to clarify the 
principle of ghost imaging, which is compatible with the 
generally used intensity fluctuation correlation theory, 
and can be used to explain the phenomena we realized 
in experiments. Moreover, the analysis of ghost imaging 
with a classical thermal source can be generalized to the 
case of an entangled source, which remains for the future 
work. 
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